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Abstract 


Strangeon stars, which are proposed to describe the nature of pulsar-like compact stars, have passed various 
observational tests. The maximum mass of a non-rotating strangeon star could be high, which implies that the remnants 
of binary strangeon star mergers could even be long-lived massive strangeon stars. We study rigidly rotating strangeon 
stars in the slowly rotating approximation, using the Lennard-Jones model for the equation of state. Rotation can 
significantly increase the maximum mass of strangeon stars with unchanged baryon numbers, enlarging the mass-range 
of long-lived strangeon stars. During spin-down after merger, the decrease of radius of the remnant will lead to the 
release of gravitational energy. Taking into account the efficiency of converting the gravitational energy luminosity to 
the observed X-ray luminosity, we find that the gravitational energy could provide an alternative energy source for the 
plateau emission of X-ray afterglow. The fitting results of X-ray plateau emission of some short gamma-ray bursts 
suggest that the magnetic dipole field strength of the remnants can be much smaller than that of expected when the 
plateau emission is powered only by spin-down luminosity of magnetars. 
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1. Introduction 


Pulsar-like compact stars, among which the well-known ones 
are radio pulsars,° are born in supernova explosions. Although 
abundant observational data of pulsars have been accumulated, 
their interior structure is still a controversial topic at present. The 
typical density of a pulsar is slightly larger than that of nuclear 
matter on average, so the equation of state (EoS) of pulsars 
essentially depends on the behavior of quantum chromodynamics 
(QCD) at low-energy scales, which is still a challenge for us to 
understand. Some basic problems still remain. Have quarks been 
de-confined there? Does strangeness play an essential role? Based 
on different points of view, a variety of models for pulsars have 
been speculated, such as neutron stars (NSs) and quark 
stars (QSs). 

A strangeon star (SS) model was proposed originally by Xu 
(2003) and has been studied for twenty years, based on both 
observational and theoretical points of view (see Lai et al. 
2023a and references therein). Briefly speaking, at realistic 
densities inside pulsars, the energy scale is much higher than 
the mass difference between strange (s) and up/down (u/d) 
quarks, while is additionally not high enough to justify the 
validity of perturbative QCD. The net strangeness could 


B They are strictly called “pulsar-like compact stars" instead of “pulsars” 
because some of them do not manifest themselves as radio pulsars. In most 
cases if it will not cause misunderstandings, they are still called “pulsars” for 
short. 


emerge due to the weak interaction, and u, d and s quarks 
would tend to be of the same amount in the light-flavor 
symmetry. At the same time, the non-perturbative effects of 
QCD could be significant, similar to the case of the strong 
interaction in nuclei. It is then conjectured that strangeons,’ 
which can be understood as “nucleons with strangeness” or 
"strange nucleons," could be the building blocks of dense 
matter in pulsars (Xu 2003; Xu & Guo 2017). 

Compact stars composed totally of strangeons are called 
strangeon stars (SSs). SSs are different from both NSs and 
strange quark stars (SQSs). Different from baryons inside NSs, 
strangeons are of three-flavored, and the number of constituent 
quarks (N4) of a strangeon can be larger than three. For 
instance, strangeons with N4 — 18 are similar to the so called 
quark-alpha (Michel 1988) which are be completely asym- 
metric in spin, flavor and color space. Moreover, NSs are 
gravitationally bound and have crusts composed of normal 
nuclei, while SSs are self-bound and have almost the same 
composition from the center to the surface. SSs are also 
different from SQSs, and the main difference is that SSs are 
self-bound by the residual interaction between strangeons, 
while SQSs are self-bound by bag-like confinement. 


Strangeons were previously called “strange quark-clusters” in some earlier 


papers. The strong coupling between quarks may group quarks into clusters 
which are composed of several quarks. 
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The SS model has passed various observational tests. It 
predicted high-mass pulsars (possibly even larger than 3M.) (Lai 
& Xu 2009a, 2009b) before the formal discovery of massive 
pulsars with M > 2M... In addition, the strangeon matter surface 
could naturally explain the pulsar magnetospheric activity (Xu 
et al. 1999) as well as the subpulse-drifting of radio pulsars (Lu 
et al. 2019). Pulsar glitch (Zhou et al. 2004, 2014; Lai et al. 
2018b) and glitch recovery (Lai et al. 2023b) can also be 
explained under the framework of starquake of SSs. The plasma 
atmosphere of SSs can reproduce the Optical/UV excess 
observed in X-ray dim isolated neutron stars(Kaplan et al. 
2011; Wang et al. 2017). The tidal deformability of merging 
binary SSs is consistent with the results of gravitational wave 
event GW170817 (Lai et al. 2018a, 2019). 

Rotation will affect the structure of pulsars, and related 
astrophysical consequences are worth exploring to provide tests 
for EoS models. A perturbative approach describing distorted 
NSs with uniform and slow rotation to the second order of 
angular frequency 22 was given by Hartle (1967) and Hartle & 
Thorne (1968), and was developed to the third order of Q by 
Hartle (1973) for calculating the variations of moments of 
inertia. Slow rotation means that Q of a star with mass M and 
radius R is much smaller than a critical frequency, 
Q « OQ, = JGM/R?, in which case the rotating configuration 
can be considered as a perturbation on a non-rotating one of the 
same central density. It has been shown that this perturbative 
approach can be applied with great accuracy for most observed 
NSs, even for most millisecond pulsars (Berti & Stergiou- 
las 2004; Benhar et al. 2005; Berti et al. 2005). Gao et al. 
(2022) have provided detailed calculations about the structure 
of slowly rotating SSs, using the EoS of the Lennard-Jones 
model, and derived the moments of inertia, the quadrupole 
moments, the eccentricities, changes in the gravitational and 
baryonic masses, and universal relations between some of these 
quantities. 

The effect of rotation on the stability is crucial for the fate of 
the products of binary NS/SS mergers. The maximum mass of 
non-rotational NSs/SSs, denoted by Mzoy, can be derived by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equations 
(Oppenheimer & Volkoff 1939) for a given EoS. According to 
the widely used definition, the rotating NS/SS that is stablized 
by differential rotation is hypermassive, the one that is stablized 
by rigid rotation is supramassive, and the one that is stable 
without rotation is stable or long-lived. The maximum mass of 
rotating NSs/SSs, denoted by Mnax, will increase with the 
angular frequency 92, which also depends on EoSs. Considering 
that the EoSs should satisfy the constraints of both the 
existence of two-solar-mass pulsars and the tidal deformability 
of GW170817 (Abbott et al. 2017), we use the EoS of the 
Lennard-Jones model (Lai & Xu 2009b) (which has also been 
used by Gao et al. 2022) for SSs and the EoS of the AP4 
model (Akmal & Pandharipande 1997) for NSs. 
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The fate of the product of a merger event should be 
determined in a constant baryon number. We will explicitly 
show the increases of Mmax with Q along the lines of constant 
baryonic mass, with the results indicating that the increases of 
Mmax are more pronounced for SSs than those for NSs. 
Combining with the previous conclusion that Mroy of SSs can 
be larger than 2.5Mo in a wide range of parameter space (Lai & 
Xu 2009a, 2009b; Lai et al. 2013; Guo et al. 2014), we can 
infer that if pulsar-like compact stars are actually SSs, the 
remnants of binary strangeon star mergers are very likely to be 
long-lived massive SSs. The long-lived SSs as the remnants of 
binary strangeon star mergers could have interesting observa- 
tional consequences, e.g., they could reproduce the light curves 
observed in kilonova(Lai et al. 2018a, 2021). As will be 
demonstrated in this paper, the long-lived massive SSs could 
also provide large gravitational energy enough to explain the 
X-ray afterglow of short gamma-ray bursts (GRBs). 

Short gamma-ray bursts (SGRBs) are generally believed to 
originate from the binary NS mergers (Eichler et al. 1989) or 
the binary NS-BH (black hole) mergers (Paczynski 1991). 
Among them, the ones with afterglow phase are often 
interpreted as originating from binary NS mergers, where the 
afterglow emission is widely accepted as being powered by the 
millisecond magnetars(e.g., Dai & Lu 1998; Zhang & 
Mészáros 2001). The electromagnetic dipolar radiation of the 
postmerger magnetars could explain the X-ray flares following 
SGRBs (Dai et al. 2006; Gao & Fan 2006), and the lightcurves 
produced by magnetar spindown winds could explain the X-ray 
plateaus at observed luminosities of SGRBs(Strang & 
Melatos 2019; Strang et al. 2021). The gravitational bursts 
due to magnetar wind dissipation of the millisecond magnetars 
left behind binary neutron star mergers (Zhang 2013) and the 
associated multi-wavelength afterglows (Gao et al. 2013) have 
been investigated. Moreover, some plateaus with long dura- 
tions are suggested to be powered by the nascent SQSs (Yu 
et al. 2009), which could be supported by observations on the 
break time of internal X-ray plateaus in SGRBs (Li et al. 2016). 
The X-ray plateaus in the afterglow of GRBs observed by Swift 
satellite have been explained by the magnetar central engines, 
where the dipole field with strength larger than 10 G is 
required (Rowlinson et al. 2013; Stratta et al. 2018). In this 
paper, we will investigate the gravitational energy released by 
long-lived massive SSs and the implication on the X-ray 
afterglow of SGRBs. 

If pulsar-like compact stars are actually SSs instead of NSs, 
i.e., the binary neutron star mergers are actually binary SS 
mergers, then SGRBs with the plateaus in the X-ray afterglow 
phase originate from binary SS mergers. It is worth noting that, 
the strong magnetic field may not be necessary for the scenario 
of SSs, e.g., the elastic and gravitational energy released of SSs 
can explain the AXPs/SGRs (anomalous X-ray pulsars/soft 
gamma repeaters) associated with glitches (Zhou et al. 2014) 
and the precursor emission of SGRBs (Zhou et al. 2023). The 
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latent heat released in the solidification of the strangeon star 
was proposed as the energy injection into the X-ray plateau of 
GRB afterglow (Dai et al. 2011), and this idea is supported by 
the X-ray light curve of GRB 170 714A whose two plateaus 
can be interpreted as being powered respectively by the latent 
heat and the spin-down of a massive strangeon star which is the 
remnant of the binary star merger (Hou et al. 2018). For the 
sake of simplicity, it was assumed that the latent heat released 
in solidification was released as the blackbody radiation to 
power the relativistic jet of GRBs; however, more realistic and 
detailed mechanism should be taken into account. 

If pulsar are actually SSs instead of NSs, it is worth 
exploring explicitly how the plateaus in the X-ray afterglows of 
SGRBs originate from binary SS mergers. In this paper, we will 
consider the contribution of gravitational energy to the X-ray 
light curves of GRBs from binary star mergers. The remnants 
of binary SS mergers will undergo spin-down due to energy 
loss. The gravitational mass M and radius R will reduce as the 
angular velocity Q reduces under the same baryonic mass My, 
so the gravitational energy will be released as an isolated star is 
spinning down. Because the remnants of binary SS mergers 
would be long-lived massive SSs, we will study the spin-down 
process of massive SSs and investigate whether the gravita- 
tional energy could provide the energy injection to the X-ray 
plateau in the afterglow of SGRBs. As we will demonstrate in 
this paper, the shrinkage of the star would lead to oscillations 
and turbulence, which would convert the gravitational energy 
into kinematic energy and finally injected into the GRB fireball. 

Assuming that the spin-down is due to magnetic dipolar 
radiation, we can derive the luminosity of gravitational energy 
releasing, then the X-ray luminosity can be derived by taking 
into account the efficiency of converting the gravitational 
energy to the observed X-ray luminosity. From Stratta et al. 
(2018) which interpreted the GRBs presenting X-Ray afterglow 
plateaus as of the magnetar origin, we choose some SGRBs 
with obvious plateaus in Swift GRB sample, and fit their X-ray 
afterglow data using the MCMC (Markov Chain Monte Carlo) 
method. Our fitting results will be compared to that in Stratta 
et al. (2018). The values of x7/ d. o. f in our scenario are not 
much larger (and in some cases are even smaller) than that in 
Stratta et al. (2018). In addition, the results show that the 
magnetic dipole field strength of the remnants can be much 
smaller than that of expected when the plateau emission is 
powered only by spin-down luminosity of magnetars. 

This paper is organized as follows. After the introduction 
about the Lennard-Jones model of SSs in Section 2.1, we 
demonstrate how to calculate the structure of slowly rotating 
SSs in the Hartle- Thorne approximation in Section 2.2, and 
show the results for maximum mass and spherical stretching 
under rotation. Based on the luminosity of gravitational energy 
during spin-down derived in Section 3.1, we investigate in 
Section 3.2 whether the gravitational energy release during 
spin-down could provide enough energy injection for the 
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plateau emission of X-ray afterglow of SGRBs. Conclusions 
and discussions are made in Section 4. 


2. Strangeon Stars in Slow Rotation Approximation 
2.1. Static Strangeon Stars 


We choose the Lennard-Jones model to describe the EoS of 
SSs(Lai & Xu 2009b) because it can well characterizes the 
non-relativistic nature and the strong-repulsive interactions at 
short distances, and the allowed parameter space to satisfy the 
constraints by both the existence of two-solar-mass pulsars and 
the tidal deformability of GW170817 is large (Lai et al. 2019). 
The potential between two strangeons is 


12 6 
atr) su] (2) - (2) | (1) 
r r 


where Up is the depth of the potential and ro is the range of 
interaction. The total energy density includes the densities of 
the potential energy, the lattice vibration energy and the 
baryonic mass-energy. The lattice energy density is negligible 
compared to the other two energy densities, so the total energy 
density is 


pc? = 2Up (6.27)? n? — 8.4r$ m?) + nmc?, (2) 
and the pressure is 
P = 4U0(12.4r}° n5 — 8.478 m), (3) 


where n is the number density of strangeons, m is the mass of 
each strangeon, and the simple-cubic lattice structure is 
assumed. If the number of quarks in each strangeon is Ng, then 
m œ Na: 300 MeV. The parameter ro is related to the baryon 
number density on the surface ny, where the pressure vanishes. 

As in Gao et al. (2022), we choose the EoS with parameters 
np s = 0.36 fm^, Uo — 30 MeV (denoted by LX3630), and 
N= 18, because it satisfies the constraint from the measure- 
ment of the moment of inertia of PSR J0737-3039A (Hu et al. 
2020). Given the EoS, the structure of non-rotating SSs is 
governed by the TOV equations (Oppenheimer & Volk- 


off 1939). The gravitational mass is M = 47 J, á pr?dr and 


" " R " 
baryonic mass is My — 4n f, py €^r)dr, where py = nm is the 


baryon density, a= 1/41 — 2Gm/(rc?) and m is the 


gravitational mass enclosed in radius r. 


2.2. Strangeon Stars under Slow Rotation 


Given a central density p, and the EoS, the structure derived 
by the TOV equation is the static and spherical background, 
based on which the gravitational mass M, the radius R at the 
equator, and the baryonic mass M, of a rigidly rotating star in 
slow rotation approximation can be derived by adding 
corrections to the second order of Q. This procedure was first 
formulated by Hartle (1967) and Hartle & Thorne (1968) for 
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Figure 1. M-R curves for SSs (red curves) and NSs (blue curves), where solid 
lines represent non-rotating configurations, and dashed lines represent rotating 
configurations with the angular frequency Q.. Sequences of constant-M, are 
denoted by dotted lines with My = MS", 24M., and 1.6Mo, where 
MEXE cs = 4.4Mo and Mf Iys = 2.7Mo. 


b,max 


rotating NSs. The structure of rotating SSs was given in detail 
by Gao et al. (2022), including the corrections induced by the 
match conditions on the surface. Here we adopt the same 
procedure as that in Gao et al. (2022) to show the evolution 
along given values of My. The values of M, M, and R are 
calculated to the spherical terms in the second order of Q, so 
only spherical deformations are considered. The moment of 
inertia J in Section 3.1 is calculated by taking into account the 
corrections to the third order of Q to the angular moment J. 

The calculation can proceed as follows, whose details can be 
found in Gao et al. (2022). (i) Choose a central density p, to 
calculate the structure of a non-rotating configuration. (ii) The 
gravitational mass M and the radius R of a rotating SS with the 
angular frequency Q. = J.GM/R? is derived by adding the 
perturbations to a non-rotating configuration under the same pe, 
taking into account the matching condition at the surface. (iii) 
The values of M and R of the configuration with a different 
angular frequency Q( < Q.) and the same central density can be 
obtained by multiplying the perturbations by the rescaling 
factor (./Q,)?. Changing the values of o, gives the M-R curve. 
The M-R curve with another value of O( < Qe.) can be derived 
from the same procedure. (iv) By connecting the same value of 
My on each M-R curve we can get a constant-M, line. 


2.2.1. Evolution under Constant Baryonic Mass 


An isolated star has an unchanged baryonic mass during 
spin-down. For a given EoS, the stable configuration with 


M = Mrov has the baryonic mass MS“ | and the sequences 
b,max 
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with baryonic mass My < Mgr will evolve to the stable 
configurations with the unchanged baryonic mass as they spin 
down. We plot the gravitational mass and radius curves of 
strangeon stars in Figure 1, under EoS LX3630 with Ng = 18, 
in which case Mss ~ 4.4M.. The sequences of 
constant-M, are denoted by the red dotted lines with 
My = 4.4Mo, 2.4Ms, and 1.6M.. For comparison, we also 
plot the result of NSs under EoS AP4 (Akmal & Pandhar- 
ipande 1997), in which case Ms ys c 2.7M., and the 
sequences of constant-M, are denoted by the blue dotted lines 
with My = 2.7M., 2.4Mc, and 1.6M.. Solid lines represent 
non-rotating configurations, and dashed lines represent rotating 
configurations with the critical angular frequency Qe. 

Our results show that the increase of gravitational mass due 
to rotation will be larger for larger mass, and the increases of 
Mmax by rigid rotation for SSs are more pronounced than that 
for NSs. The Mmax value for SSs is roughly 9.3% higher than 
Mroy along the constant baryonic mass lines, and for the case 
of NSs the result is roughly 5.6%. If only compared with the 
maximum values of the M-R curves, both the results are 
roughly 20% for SSs and NSs. However, although the increases 
in Mnax by rigid rotation of SSs are larger than that of NSs, at 
values of M well below Moy the advantage of SSs over NSs 
regarding the increased of gravitational mass M due to rotation 
over is not significant. 

In fact, the gravitational energy released by an NS during 
spin-down is usually larger than that of an SS, since the 
shrinkage of the NS is larger. Although this conclusion seems 
to favor NSs for providing gravitational energy to explain 
plateau emission, we still prefer SSs to NSs. The reason is that 
an NS as the remnant of binary neutron star mergers would not 
be long-lived. Although the AP4 model for NSs could pass 
both the tests of the massive pulsars and tidal deformability 
during mergers (Annala et al. 2018), Mroy ^ 2.2M., which is 
well below the total mass of a known binary neutron star 
system inferred by Antoniadis et al. (2016). 


2.2.2. Spherical Stretching Due to Rotation 


From Figure | we can see that between Q = Q, and Q = 0, 
the change in radius of an NS will be more than that of an SS. 
The spherical stretching due to rotation has been discussed in 
Gao et al. (2022), which was found to be increasing toward the 
surface of the star for the case of NSs while being nearly 
unchanged through the star in the case of SSs. To see the 
change of density due to rotation, we show explicitly here the 
change of central baryon density. We plot the curves showing 
the change of central baryon density p.p with Q under some 
given values of initial mass Mo = M(Q = Qo) in Figure 2, with 
the initial angular frequency Qo = 27/(1 ms). The red solid and 
dashed curves represent the results of SSs for My = 1.6M. 
(corresponding to the initial mass Mọ% 1.36M.&) and 
My — 27M. (corresponding to Mo — 2.16M.), respectively. 
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Figure 2. The change of central baryon density po» with Q for M, = 1.6M.; 
(solid) and My = 2.7M. (dashed), where red lines represent results of SSs and 


blue lines represent results of NSs, respectively. pọ is the saturated nuclear 
matter density. 


For comparison, the results of NSs are shown by the blue solid 
and dashed curves for My,— 1.6M. (corresponding to 
Mo œ 1.44Mo) and M,=2.7M. (corresponding to Mo 
2.22 M), respectively. 

From the change of ps» with Q, we can see that the spherical 
stretching of an NS is more significant than that of an SS with 
the same baryonic mass and the same initial spin frequency, 
especially when the initial mass Mo is larger than 2M... 
Conversely, during spin-down, the shrinkage of an NS is more 
significant than that of an SS, especially in the case of binary 
merger remnants. This may imply that a neutron star with 
Mo ~ 2.2Mo as the remnant of binary NS merger, although 
being long-lived, would undergo phase transition during spin- 
down. It would be interesting to explore the implication of such 
phase transition of massive or supramassive NSs, e.g., the 
energy released from the phase transition has been taken as the 
energy source related to GRBs' observations (Dai & Lu 1998; 
Sarin 2021). An SS, however, is close to the incompressible 
matter and would not undergo phase transition during 
spin-down. 

It is worth noting that, if the AP4 model applies to NSs both 
before and after mergers, the increase in central density during 
merger is much more significant than that during spin-down. 
The central density increases to almost twice as dense as that 
before merger, but only increase about 10% during spin-down. 
Therefore, if an NS with Mo ~ 2.2M. as the remnant of the 
binary NS merger would undergo phase transition during spin- 
down, the phase transition would be more likely to occur 
during merger instead of during spin-down. 
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3. Gravitational Energy Release of the Binary Merger 
Remnants 


Because Mroy values are high for strangeon stars, the 
remnants of the binary strangeon star mergers would probably 
not immediately collapse into black holes, and would be even 
long-lived. Being spinning fast at the beginning, the remaining 
massive strangeon stars will undergo spin-down due to energy 
loss. As shown in Section 2.2, the radius of a strangeon star 
will decrease as its angular frequency (2 decreases, so 
gravitational energy will be released during spin-down. It will 
be shown that, although only a fraction of gravitational energy 
will be converted into the X-ray emissions, it may play an 
important role in the X-ray afterglow of SGRBs. 


3.1. Luminosity of Gravitational Energy 


The differential rotation may be a short-term process during 
the early stages of the merger remnants, and would be no 
longer important on the longer timescales for the afterglow. In 
addition, we take the initial angular frequency Q = 27/(1 ms), 
which satisfies the slow rotation condition Q &« Q. — 
4GM/R?. Therefore, we can use the slow rotation approx- 
imation in the case of rigid rotation to derive the change of 
gravitational energy with time. 

The gravitational energy of a relativistic star can be derived 
by 


= 2 2 
Egrav = Mc* — Mpc = Exin 


R 
=4n f p — er?dr — 5/9. (4) 
0 


where the proper mass Mp is defined as Mp = 4r J, į pe^r?dr, 
and the kinetic spin energy E,j, is related to the angular 
momentum J by Ej, = JQ/2. The change of Eray with Q for a 
given My can be derived by the similar procedure used in 
Section 2.2, where the proper mass Mp is calculated to the 
second order of Q by the procedure similar to that for deriving 
M». Using the slow rotation approximation to calculate Eķin, the 
angular momentum J is calculated to the first order of Q by 
considering the rotational dragging of inertial frames, so Exin is 
also calculated to the second order of 22. 
The luminosity of gravitational energy can be derived by 


dE grav dQ 
dQ dt’ 


Leroy = Egrav = (5) 
Assuming that the spin-down is due to electromagnetic (EM) 
dipolar radiation and gravitational wave (GW) radiation, the 
change of Q with time ¢ is (Shapiro & Teukolsky 1983) 

dQ B,R«P — 32gre?05 

B : = (6) 

dt 6Ic* 5c 
where B, is the dipolar field strength at the poles, Z is the 
moment of inertia, R and c are the radius and ellipticity of the 
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star respectively. Combining Equations (5) and (6), we can get 
the luminosity of gravitational energy Lgrav. It is worth noting 
that, in the calculation of both Leray and the luminosity of 
magnetic dipole radiation Lem = B; ROA / (6c), the changes of 
R and J with €? are taken into account. R is calculated by adding 
the spherical deformation to the second order of Q, as 
demonstrated in Section 2.2. The corrections to J=J/Q is 
also to the second order of (2, since the corrections to the 
angular moment J for calculating J is derived to the third order 
of 2 (Hartle 1973; Gao et al. 2022). 


3.2. The Role of Ugray in the X-Ray Afterglow of SGRBs 


We investigate whether their gravitational energy release 
during spin-down could provide enough energy injection for 
the afterglow of SGRBs. In a supernova explosion, most of the 
gravitational energy will be taken away by neutrinos which are 
produced in phase transitions involving the weak interaction. 
However, as shown in Section 2.2.2, the shrinkage of an SS 
during the spin-down would not be large enough to cause a 
phase transition, so the loss of energy due to neutrinos would 
be unimportant. Then how will the gravitational energy 
released due to the shrinkage of an SS be injected into the 
GRB fireball? 

Similar to the process of heating the solar corona and 
accelerating the solar wind by high-frequency Alfvén 
waves (Tu & Marsch 1997; Kaghashvili 1999), the gravita- 
tional energy could be converted into kinematic one, and finally 
injected into the GRB fireball also by Alfvén waves. The 
oscillations and turbulence due to shrinkage of the star would 
lead to the magnetic reconnection and then generate Alfvén 
waves to take away the kinetic energy (which comes from 
gravitational energy), like the oscillation-driven magneto- 
spheric activity in pulsars (Lin et al. 2015). In this way, the 
form of converting gravitational energy into the fireball may 
also involve processes similar to the kinetic-energy-dominated 
shell (Zhang & Mészáros 2002) or a Poynting-flux-dominated 
outflow (Mészáros & Rees 1997). Because the efficiency ng for 
the gravitational energy to be converted into the X-ray 
emissions of GRB afterglow is unknown, all the complexities 
will be put into ng. 

In the model which interprets the X-ray afterglow plateau 
emission in GRBs as being powered by the electromagnetic 
dipolar emission from millisecond magnetars, the efficiency Nem 
of converting the dipole spin-down luminosity to the observed 
luminosity should be considered. However, jem is derived by 
fitting observational data because its prior value is hard to 
calculate, and the fitting results are usually different. Although 
some simulations suggest that Nem = 0.4—-1 (Gao et al. 2016), a 
more detailed study indicates that the X-ray radiation efficiency 
strongly depends on the saturation Lorentz factor and the typical 
value is of order 10? (Xiao & Dai 2019). Similarly, the 
efficiency ng of converting the gravitational energy luminosity to 
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the observed X-ray luminosity should also be considered, but 
how to determine 7, is still a problem. The data points of each 
SGRB that we choose are not enough to give good fitting with 
more than three free parameters, we can only choose either ng or 
Nem aS a free parameter. We find that the fitting results of setting 
Ng = Nem Will not significantly differ from that of setting jem = 0, 
so to highlight the role of gravitational energy we set nem = 0, 
except in two cases in which we set Ng = Nem to avoid rj, > 1, as 
will be shown latter. In addition, because the spin-down due to 
GW radiation would not be important in the afterglow 
phase (Zhang & Mészáros 2001), we neglect the second term of 
Equation (6). 

To test the validity of our scenario, we choose some SGRBs 
from Stratta et al. (2018) which fits a sample of GRB X-ray 
afterglows by assuming that the plateau emission is powered by 
the spin-down luminosity of millisecond magnetars. They 
derived values of y*/d. o. f of the fitting results which can be 
used to be compared with our results. Because we do not use 
model or quantitative selection criteria to identify the plateau 
phase, we choose the ones that have obvious plateaus, and 
identify the flat part of the data to be the plateau phase to begin 
our fitting. Among the ten SGRBs fitted in Stratta et al. (2018) 
we choose six ones, 051221A, 060614, 061201, 070714B, 
070809, and 090510, which have redshifts and obvious 
plateaus. Beside them, we choose another two SGRBs from 
the Swift GRB sample obvious plateaus, 130603B and 
140903A. 

For a given burst in Swift data (Evans et al. 2007, 2009), the 
source rest frame luminosity is derived from the flux F(t) at the 
time t by L(t) = 42D, (Z) F (t), where D,(z) is the luminosity 
distance at redshift z. For the correction from 1 to 10* keV in 
the burst rest frame to the observed band, the X-ray luminosity 
derived in our model is divided by the factor k,. The redshifts 
of the eight SGRBs are from Table 1 of Kisaka et al. (2017) 
(and references therein). We adopt the values of cosmological 
parameters used in Komatsu et al. (2009) to get D; (z), and ke is 
derived by the method used in Bloom et al. (2001). 

In the MCMC fits of different SGRBs, we assume 
My c 3.1M., which corresponds to Mp = 2.36M., under the 
fitted values of the initial spin period Po. Certainly, the mass 
range for remnants in binary star mergers is unknown. The total 
mass of the binary system associated with GW170817 which is 
probably larger than 2.7M (Abbott et al. 2017) and the masses 
of the known binary neutron star systems (Antoniadis et al. 
2016) can infer that, the initial mass Mo = 2.36Mo for remnants 
in binary star mergers could be reasonable. We find that the 
change of M, from 2.2Mo to Mo — 2.5M., would not make 
significant effects. 

The afterglow component from the interaction between the 
jet and interstellar medium should also be considered. We use 
the Python package “afterglowpy” (Ryan et al. 2020) which 
utilizes semianalytic approximations to the jet evolution and 
synchrotron emission to calculate afterglow light curves with 
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Figure 3. X-ray afterglow light curves of eight SGRBs with our best-fit results, fitted to the black points which begin from the flat part of the data. The afterglow 
component from the interaction between the jet and interstellar medium is shown by dashed lines. 
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Table 1 
Best-fit Results of Ours (Left) and Stratta et al. (2018) (Right) 
SGRB B, (10? G) Po (ms) Ne x7 /d. o. f. B, (10? G) x /d. o. f. 
051221A 1.697012 5.2962 0.697012 1.1017 10.4 + 0.9 1.3213 
060614 0.887013 3.33030 0.114304 1.9674 15.6 + 0.5 1.0821 
061201 4.173036 2.607031 0.14100 1.2937 20.6 + 1.5 1.2417 
070714B 8.735334 3.231050 0.584913 2.1809 36.9 + 4.0 1.9500 
070809 4.98*050 10.94*132 0.63*031 1.1828 58-11 1.2615 
090510 6.407096 2.05*025 0.25* 027 1.9234 11.6 € 0.5 1.1914 


structured jets, taking into account the relativistic beaming 
effect. In our calculations, the Gaussian jet model is used to 
calculate the contribution of interaction between the jet and 
interstellar medium, assuming that fractions of post-shock 
energy in radiating electrons €, — 0.04 and magnetic fields 
eg—-10 * (Ryan et al. 2020), the jet half-opening angle 0.. 
6.87? (Fong et al. 2015), and the truncation angle 6,, = 50.. The 
values of the number density of interstellar medium no and the 
electron power-law distribution index p are chosen from Fong 
et al. (2015) for each SGRB. The isotropic equivalent energy of 
the blast wave Ey— 1052 erg (Cao et al. 2023), and the 
viewing angle 0, can be derived from the photon index and 
peak flux for the Gaussian jet. After obtaining the flux at 3 keV 
under the above parameters using "afterglowpy," the flux at 
0.3-10 keV can be derived by the method of Gehrels et al. 
(2008). The results are shown by dashed lines in Figure 3, 
which indicate that compared with the internal energy injection, 
the afterglow component from the interaction between the jet 
and interstellar medium would not be important. 

We use the MCMC (Markov Chain Monte Carlo) method to 
fit three parameters including B,, Po, and ng. The comparison of 
the fitting results of ours and that of Stratta et al. (2018) for six 
SGRBs are shown in Table 1. The X-ray light curves 
(1-10* keV) of eight SGRBs and the histograms and contours 
for MCMC fits are shown in Figures 3 and 4, respectively. 
Because we find that the fitting results of setting Ng = Nem will 
not differ significantly from that of setting Nem = 0, we choose 
to fix Nem to avoid introducing one more parameter. To 
highlight the role of gravitational energy we set rj, — O for 
051221A, 060614, 061201, 070714B, and 070809. For 090510 
and 1409034, we set ng = Nem to avoid Ng > 1. From the fitting 
results, we can see that, for the six SGRBs from Stratta et al. 
(2018), the x /d. o. f values of our scenario are not much larger 
than that of Stratta et al. (2018), indicating that our scenario is 
at least not bad. In addition, the values of B, are much smaller 
than those expected in Stratta et al. (2018). 


4. Conclusions and Discussions 


The effects of rotation on the structure of pulsar-like 
compact stars would have interesting consequences and 
provide tests for the EoS models. Under the hypothesis that 


pulsar-like compacts stars are SSs, whose EoS is described by 
the Lennard-Jones model, we study rigidly rotating SSs in the 
slow rotation approximation. Although we only choose the EoS 
of LX3630 and N, = 18 to perform our calculations, the other 
forms of EoS under different parameters may not bring 
qualitative differences. We find that rotation can significantly 
increase the maximum mass of a stable SS. The SS with mass 
M larger than Myoy by approximately 9% can still be stable or 
long-lived during spin-down with unchanged baryonic mass 
My. Considering that Mroy of SSs can be much larger than 
2.5Mo in a wide range of parameter space, it is very likely that 
the remnants of binary SS mergers are long-lived massive SSs. 

To explore the consequences of rotating massive SSs, we 
investigate whether their gravitational energy release during 
spin-down could provide enough energy injection for the 
afterglow of SGRBs. We derive the luminosity of gravitational 
energy releasing, where the spin-down is due to magnetic 
dipolar radiation and the changes of radius R and the moment 
of inertial J with angular frequency Q are taken into account. 
The X-ray light curves can be derived by assuming that a 
fraction of gravitational energy releasing contributes to X-ray 
luminosity. By fitting X-ray afterglow of six SGRBs in Stratta 
et al. (2018) who have redshifts and obvious plateaus, we find 
that the gravitational energy released by long-lived massive 
strangeon stars could provide an alternative energy source for 
the plateau emission of X-ray afterglow. Our fitting results 
show that the magnetic dipole field strength of the remnants 
can be much smaller than that of expected in the magnetar 
scenario. The fitting results of our scenario seem not bad 
compared to the magnetar scenario which is much more 
sophisticated than ours. 

The following are some discussions. Although the plateau of 
the X-ray afterglow of SGRBs is widely accepted as being 
powered by the electromagnetic dipolar emission from 
millisecond magnetars, we demonstrate in this paper that 
gravitational energy could provide an alternative energy source. 
To avoid the complexity in the details of the millisecond 
magnetar origin scenario, we choose to fix Nem to be either 0 or 
equal to Nem. It is expected to find a reasonable way to combine 
both of them to account for the X-ray afterglow emission of 
SGRBs, by fitting a larger sample of SGRBs. 
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Figure 4. The histograms and contours for MCMC fits of the X-ray afterglow light curves of eight SGRBs in Figure 3. 
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Certainly, we use some simplifications and assumptions to 
obtain our results. The star is assumed to be rotating rigidly, 
and the slow rotation approximation is used to calculate its 
structure. Moreover, how the gravitational energy can be 
injected into GRB fireball to power the afterglow emission is 
unknown. Here we assume that a fraction of gravitational 
energy could be converted into kinematic one and finally 
injected into the GRB fireball by Alfvén waves, which could be 
generated by oscillations and turbulence due to the shrinkage of 
the star, and we put all the ignorance into the efficiency ng. 
Similar to the efficiency nem in the magnetar scenario, ng would 
not be a constant and would depend on many factors, such as 
the injected luminosity and the process of injection. The 
improved version of our scenario in the future by fitting a larger 
sample of SGRBs would be promising. 

Some further investigations about GRB afterglow are 
expected. We only consider in this paper the afterglow of the 
SGRBs instead of the long GRBs (LGRBs), which are 
generally believed to originate from the supernova explosions 
whose remnants are pulsar-like compact stars or black holes. 
The millisecond magnetars, if the magnetic dipole field could 
be high enough, are generally accepted to be the engine of 
X-ray afterglow plateaus for both LGRBs and SGRBs. 
However, for the case of LGRBs whose remnants are stable 
compact stars with mass of about 1.4Mo, the gravitational 
energy released during spin-down alone may not be large 
enough to account for the X-ray afterglow plateaus, no matter 
whether the remnants are SSs or NSs. The real process of 
energy injection into the GRB afterglow would be complex, 
and the gravitational energy would only be a part of the 
available energy sources. 

In addition, we only consider the long-lived remnants of 
binary mergers and their implications on the X-ray afterglow 
plateaus of SGRBs. Some X-ray afterglows of SGRBs show 
rapid decay after the plateau phase, in which case the plateau is 
interpreted to be powered by the supramassive remnants and 
the rapid decay is thought to be the signal of collapse into BHs. 
The study about the distribution of the break time (i.e., collapse 
time) of X-ray plateau in SGRBs could provide test for EoS 
models of NSs and QSs (Li et al. 2016). For SS model, the 
implication for the break time of X-ray plateaus in SGRBs is 
worth exploring. On one hand, to calculate the break time we 
should know the mass-distribution of the binary neutron stars 
via population synthesis and more detailed observations on 
SGRBs, which at present have many uncertainties. On the other 
hand, the sudden decrease of X-ray light curves may not imply 
the collapse into black holes, since the decay of L4, could be 
rapid even for a long-live massive SS. As the first attempt to 
explore the related issues, in this paper, we only consider the 
case of long-lived remnants and try to explain the observed 
luminosity of X-ray plateaus. Further works about how to test 
different EoS models by observations of GRBs are expected. 


10 


Yang et al. 


Acknowledgments 


We would like to thank Dr. Yong Gao and Prof. Yun-Wei 
Yu for useful discussions, and to thank the anonymous referee 
for comments. This work makes use of data supplied by the UK 
Swift Science Data Centre at the University of Leicester. This 
work is supported by the National SKA Program of China 
(Nos. 2020SKA0120300, 2020SKA0120100), the Outstanding 
Young and Middle-aged Science and Technology Innovation 
Teams from Hubei colleges and universities (No. T2021026), 
the Young Top-notch Talent Cultivation Program of Hubei 
Province, and the Key Laboratory Opening Fund (MOE) of 
China (grant No. QLPL2021P01). 


References 


Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, PhRvL, 119, 161101 

Akmal, A., & Pandharipande, V. R. 1997, PhRvC, 56, 2261 

Annala, E., Gorda, T., Kurkela, A., & Vuorinen, A. 2018, PhRvL, 120, 172703 

Antoniadis, J., Tauris, T. M., Ozel, F., et al. 2016, arXiv:1605.01665 

Benhar, O., Ferrari, V., Gualtieri, L., & Marassi, S. 2005, PhRvD, 72, 044028 

Berti, E., & Stergioulas, N. 2004, MNRAS, 350, 1416 

Berti, E., White, F., Maniopoulou, A., & Bruni, M. 2005, MNRAS, 358, 923 

Bloom, J. S., Frail, D. A., & Sari, R. 2001, AJ, 121, 2879 

Cao, X.-F., Tan, W.-W., Yu, Y.-W., & Zhang, Z.-D. 2023, arXiv:2306.16795 

Dai, S., Li, L., & Xu, R. 2011, SCPMA, 54, 1541 

Dai, Z. G., & Lu, T. 1998, A&A, 333, L87 

Dai, Z. G., Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 311, 1127 

Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Natur, 340, 126 

Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 379 

Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177 

Fong, W., Berger, E., Margutti, R., & Zauderer, B. A. 2015, ApJ, 815, 102 

Gao, H., Ding, X., Wu, X.-F., Zhang, B., & Dai, Z.-G. 2013, ApJ, 771, 86 

Gao, H., Zhang, B., & Lü, H.-J. 2016, PhRvD, 93, 044065 

Gao, W.-H., & Fan, Y.-Z. 2006, ChJAA, 6, 513 

Gao, Y., Lai, X.-Y., Shao, L., & Xu, R.-X. 2022, MNRAS, 509, 2758 

Gehrels, N., Barthelmy, S. D., Burrows, D. N., et al. 2008, ApJ, 689, 1161 

Guo, Y.-J., Lai, X.-Y., & Xu, R.-X. 2014, ChPhC, 38, 055101 

Hartle, J. B. 1967, ApJ, 150, 1005 

Hartle, J. B. 1973, Astrophysics and Space Science, 24, 385 

Hartle, J. B., & Thorne, K. S. 1968, ApJ, 153, 807 

Hou, S.-J., Liu, T., Xu, R.-X., et al. 2018, ApJ, 854, 104 

Hu, H., Kramer, M., Wex, N., Champion, D. J., & Kehl, M. S. 2020, MNRAS, 
497, 3118 

Kaghashvili, E. K. 1999, Geophys. Res. Lett., 26, 1817 

Kaplan, D. L., Kamble, A., van Kerkwijk, M. H., & Ho, W. C. G. 2011, ApJ, 
136, 117 

Kisaka, S., Ioka, K., & Sakamoto, T. 2017, ApJ, 846, 142 

Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 

Lai, X., Xia, C., & Xu, R. 2023a, AdPhX, 8, 2137433 

Lai, X.-Y., Gao, C.-Y., & Xu, R.-X. 2013, MNRAS, 431, 3282 

Lai, X. Y., Wang, W. H., Yuan, J. P., et al. 2023b, MNRAS, 523, 3967 

Lai, X.-Y., Xia, C.-J., Yu, Y.-W., & Xu, R.-X. 2021, RAA, 21, 250 

Lai, X.-Y., & Xu, R.-X. 2009a, APh, 31, 128 

Lai, X.-Y., & Xu, R.-X. 2009b, MNRAS, 398, L31 

Lai, X.-Y., Yu, Y.-W., Zhou, E.-P., Li, Y.-Y., & Xu, R.-X. 2018a, RAA, 
18, 024 

Lai, X.-Y., Yun, C.-A., Lu, J.-G., et al. 2018b, MNRAS, 476, 3303 

Lai, X.-Y., Zhou, E.-P., & Xu, R.-X. 2019, EPJA, 55, 60 

Li, A., Zhang, B., Zhang, N.-B., et al. 2016, PhRvD, 94, 083010 

Lin, M.-X., Xu, R.-X., & Zhang, B. 2015, ApJ, 799, 152 

Lu, J.-G., Peng, B., Xu, R.-X., et al. 2019, SCPMA, 62, 959505 

Mészaros, P., & Rees, M. J. 1997, ApJL, 482, L29 

Michel, F. 1988, PhRvL, 60, 677 

Oppenheimer, J. R., & Volkoff, G. M. 1939, PhRv, 55, 374 

Paczynski, B. 1991, AcA, 41, 157 


Research in Astronomy and Astrophysics, 24:035005 (11pp), 2024 March 


Rowlinson, A., O'Brien, P. T., Metzger, B. D., Tanvir, N. R., & Levan, A. J. 
2013, MNRAS, 430, 1061 

Ryan, G., van Eerten, H., Piro, L., & Troja, E. 2020, ApJ, 896, 166 

Sarin, N. 2021, PhD thesis, Monash University, Australia 

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs and 
Neutron Stars. The Physics of Compact Objects (New York: Wiley) 

Strang, L. C., & Melatos, A. 2019, MNRAS, 487, 5010 

Strang, L. C., Melatos, A., Sarin, N., & Lasky, P. D. 2021, MNRAS, 507, 2843 

Stratta, G., Dainotti, M. G., Dall'Osso, S., Hernandez, X., & De Cesare, G. 
2018, ApJ, 869, 155 

Tu, C. Y., & Marsch, E. 1997, SoPh, 171, 363 

Wang, W.-Y., Lu, J.-G., Tong, H., et al. 2017, ApJ, 837, 81 

Xiao, D., & Dai, Z.-G. 2019, ApJ, 878, 62 


11 


Yang et al. 


Xu, R., & Guo, Y. 2017, in Centennial of General Relativity: A Celebration., 
ed. C. A. Zen Vasconcellos (Singapore: World Scientific Publishing Co. 
Pte. Ltd.), 119 

Xu, R.-X. 2003, ApJL, 596, L59 

Xu, R.-X., Qiao, G.-J., & Zhang, B. 1999, ApJL, 522, L109 

Yu, Y.-W., Cao, X.-F., & Zheng, X.-P. 2009, ApJL, 706, L221 

Zhang, B. 2013, ApJL, 763, L22 

Zhang, B., & Mészáros, P. 2001, ApJL, 552, L35 

Zhang, B., & Mészáros, P. 2002, ApJ, 566, 712 

Zhou, A.-Z., Xu, R.-X., Wu, X.-J., & Wang, N. 2004, APh, 22, 73 

Zhou, E., Gao, Y., Zhou, Y., et al. 2023, arXiv:2305.10682 

Zhou, E-P. Lu, J.-G., Tong, H., & Xu, R.-X. 2014, MNRAS, 443, 
2705 


